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(54) Multi-user detection CDMA receivers in mobile telecommunications systems 



(57) In order to reduce complexity of computation of 
Multiple Access interference (MAI) in a multi-user de- 
tection receiver of a mobile telecommunications system, 
at least some of the cross- correlations (p) between user 
spreading codes (c) of the correlation matrix are com- 
puted (20) during the call set up process prior to call 
initiation, and such cross-code correlations are stored 
in a look up table (22) for use in calculation of the cor- 
relation matrix during reception of the user signals. The 
multi-user detection receiver includes means (10) for 



forming a channel estimate (h) for each of a plurality of 
users (1 ...K), means (12, 1 4) for forming a system matrix 
(A), and applying the system matrix to a rake filter (14) 
to modify incoming signals (r), means (18) for forming a 
correlation matrix (R), and multi-user detection (MUD) 
means (1 6) employing the correlation matrix to recover 
user symbols. The invention is applicable both to joint 
and to symbol-level subtractive interference cancella- 
tion detectors, and more generally to any receiver em- 
ploying multi-user parameter estimation. 
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Description 

[0001] The present invention relates to Code Division Multiple Access (CDMA) receivers in mobile telecommunica- 
tions systems, particularly though not exclusively multi-user detection receivers. 
5 [0002] The issues involved in satisfactory reception at a Base Transceiver Station (BTS) of multiple transmissions 
from mobile stations (MS), modulated with spread spectrum "spreading codes" are well documented, see the following 
references: 

[1] S. Verdu, "Multiuser Detection", Cambridge University Press, 1998, ISBN 0 521 59373 5. 
10 [2] S. Moshavi, "Multi-User Detection for DS-CDMA Communications", IEEE Communications Magazine, pp. 

124-136, October 1996. 

[3]A. Duel-Hallen, J. Holtzman, Z.Zvonar, "Multiuser Detection for CDMA Systems", IEEE Personal Communica- 
tions Magazine, pp. 46-58, April 1995. 

15 [0003] As appears from these references, the spreading codes are in principle orthogonal, and should not interfere 
with one another to enable detection by "matched filter" correlation . In practice however for a variety of reasons, e.g. 
time delay between the various received signals, there exists cross correlation between the received signals - multiple 
access interference (MAI). This greatly complicates satisfactory reception, the mathematical complexities making an 
optimal detector impractical. Various approximation techniques have therefore been devised such as decorrelating 

20 detectors which applies the inverse of the correlation matrix to the output of the matched filters. Other detectors are 
the class of linear detectors including the minimum mean-squared error (MMSE) detector, which minimises the mean 
squared error between the actual data and the soft output of the detector, and the polynomial expansion (PE) detector 
which applies a polynomial expansion in the correlation matrix to the matched filter bank output. There is also a class 
of detectors employing subtractive interference cancellation which creates at the receiver separate estimates of the 

25 MAI contributed by each user in order to subtract out some or all of the MAI seen by each user. In order to simplify 
computation, it has been proposed in [1], that knowledge of or processing of cross correlations may be avoided in 
certain circumstances, and Reference [2] points out with short spreading the partial correlations for each bit can be 
assumed the same. 

[0004] Despite the large number of detectors which have been proposed and implemented, there remains a need 
30 for improved detectors for multi-user detection CDMA receivers in mobile telecommunications systems. 

SUMMARY OF THE INVENTION 

[0005] The concept of the present invention is to provide a receiver which enables a significant part of the computation 
35 of the correlation matrix to be carried out prior to the establishment of a communications link. This therefore enables 
a reduced complexity of computation during the reception process, (or alternatively enables use of a more complex 
algorithm without increasing overall complexity of real-time computation). Thus computational complexity associated 
with the computation of the correlation matrix in multi-user parameter estimation in a CDMA system, specifically multi- 
user detection (MUD) CDMA receiver, is reduced. 
40 [0006] A time-sampled signal segment received at a base-transceiver station (BTS)can be written as a vector r = 
Ax+n where r is the vector of unknown symbols transmitted by all K mobile users in the cell, n is the noise vector. Matrix 
A is a sparse multi-user "system matrix" whose columns represent the channel-filtered user spreading codes. The 
objective of multiuser detection is to estimate the vector x of all user symbols, given the received vector r and the matrix 
A. The implementation of a MUD scheme according to the invention (including joint-detection and symbol-level inter- 
ns ference cancellation algorithms) involves the computation of a correlation matrix R=A H A whose elements represent 
the cross-correlations between various scaled and delayed replicas of the users' spreading codes. The computation 
of R, which significantly contributes to the overall computational complexity, is reduced by noting in accordance with 
the invention that the user spreading codes are known a priori at the base station. As a result, the code cross-corre- 
lations are computed off-line and stored prior to the establishment of the communication link. The stored values are 
so retrieved in real time to implement the multi-user detection process. 

[0007] Thus in a first aspect, the invention provides a Code Division Multiple Access (CDMA) receiver for a mobile 
telecommunications system, including means for forming a channel estimate for each of a plurality of users, means 
for forming a system matrix wherein rows and columns represent users and spreading codes as modified by channel 
estimates , means for forming a correlation matrix, and means employing said correlation matrix and system matrix to 
55 recover user symbols, characterised in that said means for forming a correlation matrix includes means for computing 
prior to signal reception at least some cross-correlations between user spreading codes , and storage means for storing 
such cross-code correlations for use in calculation of the correlation matrix. 

[0008] The invention is applicable in principle to any type of CDMA receiver having a multi-user detector employing 
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a correlation matrix, specifically joint detectors and symbol-level subtractlve interference cancellation detectors. 
[0009] More generally, the invention relates to reduction of computational complexity of multi-user parameter esti- 
mation in CDMA systems, and accordingly the invention provides in a more general aspect, a Code Division Multiple 
Access receiver for a mobile telecommunications system, including means for multi-user parameter estimation for a 
5 plurality of users, including means for forming a correlation matrix, characterised in that said means for forming a 
correlation matrix includes means for computing prior to signal reception at least some cross- correlations between 
user spreading codes, and storage means for storing such cross-code correlations for use in calculation of the corre- 
lation matrix. 

[0010] In a further aspect, the invention provides in a multi-user detection Code Division Multiple Access (CDMA) 
10 mobile telecommunications system, a method of receiving individual user signals, comprising: 

forming a multi-user parameter estimation for a plurality of users (1 ...K), including forming a correlation matrix (R), 
and employing said correlation matrix to recover user symbols, characterised in that the method comprises de- 
tecting a new call, recovering the spreading code for the respective kth mobile station at the base transceiver 
15 station (BTS) during the call set-up process; 

computing during the call set-up process cross-correlations between the codes of mobile station k and those of 
the k-1 mobile stations which have already established calls with the BTS; and 

employing the cross-correlations during the call for reception of the symbols transmitted by all k mobile stations. 

20 BRIEF DESCRIPTION OF THE DRAWINGS 

[0011] A preferred embodiment of the present invention will now be described with reference to the accompanying 
drawings, wherein:- 

25 Figure 1 is a diagram of a system matrix employed in the present invention; 

Figure 2 is a block schematic diagram of a known multi-user detection receiver for multi-user CDMA signals; 
Figure 3 is a block schematic diagram of the preferred embodiment of the multi-user detection receiver according 
to the invention; 

Figure 4 is a diagram of a correlation matrix employed in the present invention; 
30 Figures 5 and 6 are explanatory diagrams showing the computation of sub-matrices of the correlation matrix; and 

Figure 7 is a flow chart showing a method of call establishment in accordance with the invention. 

DESCRIPTION OF THE PREFERRED EMBODIMENT 

35 [001 2] A time-sampled segment of the signal received at a CDMA base-station can be written as a complex vector 
r = Ax + n where x is the vector of unknown complex symbols transmitted by all K mobile users in the cell and n is the 
complex noise vector. The complex matrix A is a multi-user "system matrix 11 whose columns represent the user spread- 
ing codes which are modified by the filtering action of the radio channel. 

[0013] Without loss of generality (users may have different spreading factors, i.e. different data rates), consider a 
40 scenario as shown in the matrix representation of Figure 1 involving K cochannel asynchronous mobile transmissions 
each using a spreading factor of Q. The first K columns represent K transmissions of the first symbol, and the nth K 
columns represent the nth symbol for each of K users. In each column, there is contained the user spreading code as 
modified by the respective user channel. Thus the system matrix may be represented by an array of sub-matrices B, 
one sub-matrix for each user symbol, generally arranged along the diagonal of A. Figure 1 shows the system matrix 
45 for a chip-sampled time slot duration corresponding to N bits per user. Each sub matrix B may be written as 



M 

where the non-zero elements of the k* h column, b< k >, of sub-matrix B represent the convolution between the UP* user's 
complex spreading code p< k ) (of length Q) and the k th user's propagation channel h(*) or channel estimate (of length 
55 W). As is well known in the art, the array of sub-matrices B, one sub-matrix for each user symbol, generally arranged 
along the diagonal of A. Figure 1 shows the system matrix for a chip-sampled time slot duration corresponding to N 
bits per user. Each sub matrix B may be written as where the non-zero elements of the k to column, fc>00, of sub-matrix 
B represent the convolution between the k* user's complex spreading code c< k ) (of length Q) and the k m user's prop- 
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agation channel h< k > or channel estimate (of length W). As is well known In the art, the channel hM is typically repre- 
sented by a filterhaving L taps (or fingers) of complex amplitudes tflf 1=1 ...L and delays xf^ < x ( $ < ... < x^. The 
channel is assumed fixed over a time slot period, and the notation c< k > f^W) represents the code cM delayed by W 
Many multi-user detection schemes involve the computation of the correlation matrix R=A H A. A H is the complex trans- 
5 posed skew symmetric representation of the matrix A, and the correlation matrix R is Hermitian. It will be appreciated 
that other representations of the correlation matrix are possible, and that the present invention does not depend on 
the precise formulation of the correlation matrix. 

[0014] The block diagram for a known multi-user-detection receiver, embodying the above considerations, is depicted 
in Figure 2. In figure 2, an incoming signal r is applied to a channel estimation unit 10. Typically this forms for each 

10 user from an initial training sequence the coefficients h^*), x^. These coefficients are applied to a system matrix 
computation unit 1 2, which, together with the user spreading codes cW computes the matrix A. The complex transposed 
skew symmetric form of the matrix, AH is applied to a rake filter 14 , having L taps or fingers. The incoming signal is 
modified in rake filter 1 4 according to A H . The modified signal y is applied to multi-user detection unit 1 6. The correlation 
matrix A H A is computed in unit 18 and applied to MUD unit 16, for recovering the user symbols. Thus the signal y 

is applied to MUD unit 1 6 from rake filter 1 4 includes cross correlation interference (MAI). It is the task of unit 1 6 to remove 
such cross correlations, and derive estimates of the symbols x transmitted by the users. 

[0015] In a Joint detection receiver, MUD unit 16 derives joint estimates of symbols x by simultaneously solving 
the equations for the signal y, in known manner. It will be noted the present invention is also applicable to the known 
class of symbol-level subtractlve interference cancellation detectors. These create at the receiver separate esti- 
va mates of the MAI contributed by each user in order to subtract out some or all of the MAI seen by each user. Both in 
the case of the joint detector and the subtractive interference cancellation detector, computation of the correlation 
matrix R is necessary. 

[0016] Referring now to Figure 3, the preferred multi-user detection receiver of the present invention is shown in 
block schematic form. Similar parts to those of Figure 2 are represented by the same reference numerals. In figure 3, 

25 unit 20 is provided for computing spreading code cross correlations p prior to initiation of a user call, as will be described 
more specifically below. Only spreading codes of active users are applied to unit 20, and to system matrix unit 12. 
These spreading code cross correlations p are stored in a look up table 22. The look up table 22 is accessed by the 
channel delay values x^ k ), computed by channel estimate unit 1 0, to provide the appropriate values of cross correlations 
p to the correlation matrix unit 18. Thus unit 18 is relieved of the burden of computing the cross correlation values p 

30 during a user call. It will be observed that cross correlation unit 18 requires in addition the values h-jW computed by 
estimate unit 1 0. It does not require the system matrix computed in unit 1 2, in contrast to the arrangement of Figure 2. 
[0017] In more detail, the correlation matrix of the present invention R=A H A computed in unit 1 8 may be represented 
as shown in Figure 4, wherein sub matrices B in A result in sub matrices in R, comprising matrices R ti1 and cross 
matrices R^ and Rgj" The effort required for the computation of R isas follows. 

35 

COMPUTATION OF R 1fl 

[0018] The sub-matrix l? 1f1 = B"0 is Hermitian (see above) . The computation of R 1(1 is shown diagrammatically in 
Figure 5, as multiplication of B with its transposed skew symmetric form. The resulting coefficients are as shown in the 
40 lower diagonal of the matrix, with K(K+1) unique elements given by 

kil -ft^V-fc t h<'W«c l '»(T<>>)c M (T?) 

pr r=r<'>-r?> 

where the computation of each inner product pM requires an average of VfcQ multiply-and-accumulates of the form 
(±1±j1)(±1ij1) + ... = (±1±1) + j(±1±1) + ... implying a total of %QX4 negations (NEGS) and l*Qx4 real additions 
so (RADDS). The computation of each product-and-add tyP) x Hf*? x p*<* + ... requires 1 6 real multiplications RMULTS 
and 8 real additions RADDS. It will be appreciated that NEGS, RADDS and RMULTS are software instructions each 
of a single machine cycle in a DSP. 
The total complexity for computing R A , is given by 

55 
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RADDS 
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RMULTS 




•/JC(K+1) 16L Z 



Table 1 - Operation count for computation of Ri.i. 



[0019] For implementation on a typical programmable DSP where a negation, real addition and real multiplication 
are each single-cycle instructions, the MIPS count would be as follows: 





Computation of p 


Other 


NEGS 


9.1238xl0 3 MIPS 




RADDS 


9.1238 xl0 3 MIPS 


0.2852 xlO 3 MIPS 


RMULTS 




0.5702x10* MIPS 



Table 2 - Programmable DSP MIPS count for a time slot of 667 |is 
with K-54, Q=128. 



COMPUTATION OF Rj> , 

[0020] The sub-matrix R^i equals the inner product of one segment, B 1 of B with a different segment, of B. Or 
more specifically 

fl 2,1 =B £ B 1 = [ 0 ll!a+IV-2 l 1:K[ S laf12O+IV-2 i 1:K 

as depicted in Figure 6. It will be seen the multiplication of B 1 with the skew symmetric transposed matrix B2 H results 
in a truncation of the spreading codes of both sub-matrices with a resulting square matrix Ffe j having K 2 elements. 
Thus we are now dealing with inner products between delayed and "truncated replicas" of user spreading codes. The 
extent of truncation for each spreading code is a function of the its corresponding channel delays and needs to be 
accounted for when computing the inner products. The matrix Rg ^ is not symmetric and has K 2 unique elements given 
by 

Px T ~ T I T X 

where £ represents a truncated version of v. The computation of each inner products requires an average of approx- 
imately 16x160 multiply-and-accumulates of the form (±1± j1)(±1±j1)+ ... = (±1±1) + j(±1±1) + ... implying a total of 
Hx%Qx4 NEGS and S6x%Qx4 RADDS. The computation of each product-and-add ttf) X x#* + ... requires 
1 6 RMULTS and 8 RADDS. The total complexity for computing F^ t1 is given by 
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K 2 L Z (Q) 




RADDS 


K 2 L 2 (Q) 


K 2 (8L 2 ) 


RMULTS 




K 2 (16L 2 ) 



10 Table 3 - Operation count for computation of R2 ,1 . 



For implementation on a typical programmable DSP where a negation, real addition and real multiplication are each 
15 single-cycle instructions, the MIPS count would be as follows: 





Computation of p 


Other 


NEGS 


8.9580x10* MIPS 




RADDS 


8.9580 xlO 3 MIPS 


0.5599x10* MIPS 


RMULTS 




1.1197 xlO 3 MIPS 



Table 4 - Programmable DSP MIPS count for a time slot of 667 us 
with K=54, 0=128. 



[0021] Thus, the computation of the R t 1 and R2 t i correlation sub-matrices requires the computation of the inner- 
products between delayed replicas of the user's spreading codes. This is a chip-level operation and, as is evident from 
Tables (2) and (4), contributes significantly to the overall computational burden. Since the base-station receiver is, by 
35 definition, aware of the spreading codes used by the active mobiles In its cell, it is possible, as shown in Figure 3 to 
compute the inner products off-line prior to the establishment of the communication session. Since the code delays 
are a function of radio channels which vary significantly during the course of a communication session, it is necessary 
to pre-compute the inner products for all possible combinations of multipath delays. 

[0022] Specifically, regarding R 1t , the inner products pf * p,q=1...K 1-fe=-(Q-1)...(Q-1) need to be pre-computed 

40 and stored in memory look-up table 22. For K users, there are a total of K 2 possible combinations of p and q. However, 
only %K(K+1) of these are required since R 1 f1 is symmetric. For each combination of p and q there are 2Q-1 possible 
relative code delays. The extreme values of the real and imaginary components of are +2Q and -2Q respectively. 
Then assuming O n)8 x= 1 28 , two 1 0-bit words are sufficient for storing each value of p resulting in a total memory storage 
requirement of #K(K+1) x (2Q-1) x 20 bits. For K=54 and Q=128, this implies - 1 Mbytes of memory for the look-up 

45 table. This amount of memory provides a saving by a factor of about X22 in a DSP-based computation of R 1(1 . Re- 
garding Rgj, the inner products p?* p,q=1 ...K 1 -A^=-(Q-1)...(Q-1) need to be pre-computed and stored in another region 
of memory look-up table 22. For K users, there are a total of K 2 possible combinations of p and q. For each combination 
of p and q there are 2Q-1 possible relative code delays (many of whom result in zero inner products due to the truncation, 
however this will be ignored). The extreme values of the real and imaginary components of pP* are +2Q and -2Q 

so respectively. Then assuming Qmax^l 28, two 10-bit words are sufficient for storing each value of p resulting in a total 
memory storage requirement of K 2 x (2Q-1) x 20 bits. For K=54 and Q=128, this implies — 2Mbytes of memory. This 
amount of memory provides a saving by a factor of about x 12 in a DSP-based computation of R^ . 
[0023] Note that the look-up table 22 need only contain the code cross-correlations of the mobiles that have already 
established dedicated communication channels on the uplink. 

55 [0024] A sequence of events for the use of the embodiment of Figure 3 in computation and storage of code cross- 
correlations is shown in the flow chart of Figure 7, as follows. 

[0025] In an initial step 70, k-1 mobiles stations (MS) have already established dedicated communication channels 
with a BTS (Base Transceiver Station) to provide a signal r. 
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[0026] A repetitive check is made at 72 for a new call originating from a mobile via a random-access channel or 
the BTS initiating a network-originated call via a paging channel targetted at the mobile. 
[0027] In either case, the mobile specifies the nature of its spreading code, c( k ), to the BTS during the call set-up 
process, as at 74. The BTS then computes the cross-correlations and J%* q=1 ...k-1 between the codes of mobile 
5 k and those of the k-1 mobiles which have already established dedicated chan nets with the BTS. The cross-correlations 
are stored in look up table 22. 

[0028] In step 76, during the call, relevant values are accessed from the look up table by addressing the look up 
table with the delay values x^ k > computed in unit 1 0, and the retrieved values from the look up table are applied to the 
correlation matrix unitl 8. Values of the channel estimate hjM are applied from estimator unit 10, and the sub matrices 
10 R 1 1( R 2 j of the correlation matrix are computed, as explained above. The matrix R is then applied to MUD unit 1 6 to 
perform multi-user detection of the symbols transmitted by ail k mobiles. 

[0029] A repetitive check is made at 78 to determine whether the call for user k is terminated. Once the call between 
the base and user k is determined as terminated as at 80, the cross-correlation values for user k are erased from the 
look-up table since they are no longer needed. 

15 

Claims 

1. A Code Division Multiple Access (CDMA) receiver for a mobile telecommunications system, including means 
20 (1 r>i 8) for multi-user parameter estimation for a plurality of users (1 ...K), including means for forming a correlation 

matrix (R), characterised in that said means for forming a correlation matrix includes means (20) for computing 
prior to signal reception at least some cross- correlations (p) between user spreading codes (c), and storage means 
(22) for storing such cross-code correlations for use in calculation of the correlation matrix. 

25 2. A multi-user detection Code Division Multiple Access (CDMA) receiver for a mobile telecommunications system, 
including means (10) for forming a channel estimate (h) for each of a plurality of users (1 ...K), means (12, 14) for 
forming a system matrix (A) wherein rows and columns represent users and spreading codes (c) as modified by 
channel estimates (h), means (18) for forming a correlation matrix (R), and means (14, 16) employing said corre- 
lation matrix and system matrix to recover user symbols, characterised in that said means for forming a correlation 

30 matrix includes means (20) for computing prior to signal reception at least some cross- correlations (p) between 

user spreading codes (c), and storage means (22) for storing such cross-code correlations for use in calculation 
of the correlation matrix. 

3. A receiver according to claim 2, including rake filter means (14) for applying a version of the system matrix to 
35 incoming signals (r), and mufti-user detection means (1 6) employing said correlation matrix to recover usersymbols 

4. A receiver according to any preceding claim, wherein the storage means comprises a look-up table (22), and 
wherein the cross-correlation computing means is arranged to compute all possible cross correlations for storage 
in the look up table. 

40 

5. A receiver according to claim 2, wherein the storage means comprises a look-up table (22), and wherein the cross- 
correlation computing means is arranged to compute all possible cross correlations for storage in the look up table, 
and wherein the look up table is addressed by said channel estimate means so as to provide cross correlations 
for estimated channel delays (x) to said means for forming a correlation matrix. 

45 

6. A receiver according to claim 2 or 5, wherein said means for forming a correlation matrix is arranged to receive 
channel estimate coefficients (h) from said means for forming a channel estimate, and is arranged for multiplying 
these coefficients by the cross-code correlations. 

50 7. A receiver according to any preceding claim, wherein the means for forming a correlation matrix is effective to 
divide the system matrix into submatrices (B), one for each user symbol of a message of N symbols, and wherein 
each submatrix is multiplied by a version of an adjacent sub matrix to provide a further sub matrix comprising 
diagonal sections R 1 (1 and cross sections , R^, R2 t M . 

55 8. A receiver according to claim 2, wherein the means for forming a correlation matrix is effective to divide the system 
matrix into submatrices (B), one for each user symbol of a message of N symbols, and wherein each submatrix 
is multiplied by a version of an adjacent sub matrix to provide a further sub matrix comprising diagonal sections 
R 1t and cross sections , R 21 , R^H and wherein the storage means is arranged to provide cross correlation coef- 



7 



EP 1 128 564 A1 



ficients pj*, q=1 ...k-1 , of a diagonal section for channel delay values t as follows: 



where h represents channel estimates 

9. A receiver according to claim 2 or 8, wherein the means for forming a correlation matrix is effective to divide the 
system matrix into submatrices (B), one for each user symbol of a message of N symbols, and wherein each 
submatrix is multiplied by a version of an adjacent sub matrix to provide a further sub matrix comprising diagonal 
sections R„ and cross sections , R^ , R2i H f and wherein the storage means is arranged to provide truncated cross 
correlation coefficients fit* q=1 ...k-1 , of a cross section R2 f1 for channel delay values T as follows 



where h represents channel estimates 

10. A receiver according to claim 2, wherein the receiver is a joint detection receiver 

11 . A receiver according to claim 2, wherein the receiver includes a subtractive interference cancellation detector. 

12. In a multi-user detection Code Division Multiple Access (CDMA) mobile telecommunications system, a method of 
receiving individual user signals, comprising: forming a multi-user parameter estimation for a plurality of users (1 ... 
K), including forming a con-elation matrix (R) and employing correlation matrix to recover user symbols, charac- 
terised in that the method comprises detecting a new call, recovering the spreading code( c< k >)f or the respective 
mobile station (k) at the base transceiver station (BTS) during the call set-up process; 

computing during the call set-up process cross-correlations (p^ Q , # f ) between the codes of mobile station k 

and those of the k-1 mobile stations which have already established calls with the BTS; and 

employing the cross-correlations during the call for reception of the symbols transmitted by all k mobile stations. 

13. A method according to claim 12, wherein the new call originates from a k* mobile station via a random-access 
channel or the BTS initiates a network-originated call via a paging channel targetted at the k* h mobile station. 

14. A method according to claim 12 or 13, including detecting termination of the call, and once the call is terminated, 
erasing the relevant cross correlation values from the storage means. 

15. A method according to claim 12, 13 or 14, wherein all possible cross correlations are stored in a look up table (22), 
and the look up table is addressed by estimated channel delays (x), so as to provide cross code correlations for 
forming a correlation matrix. 

16. A method according to claim 15, wherein the correlation matrix is formed from channel estimate coefficients (h), 
and multiplying these coefficients by the cross-code correlations. 

17. A method according to any of claims 12 to 16, wherein the correlation matrix is formed by dividing the system 
matrix into submatrices (B), one for each user symbol of a message of N symbols, and wherein each submatrix 
is multiplied by a version of an adjacent sub matrix to provide a further sub matrix comprising diagonal sections 
R 1t1 and cross sections , fy,^ R2,i H - 
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18. A method according to claim 17, wherein the cross correlation coefficients p*« q=t ...k-1 , of a diagonal section 
Rj t1 for channel delay values i are as follows: 

where h represents channel estimates. 

19. A method according to claim17 or 18, wherein truncated cross correlation coefficients ft* , q=1 ...k-1 , of a cross 
section for channel delay values x are as follows 

£ ± ftW £""(*!")£"'(*?') ■ 

where h represents channel estimates. 
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FIG. 5 
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FIG. 7 




Obtain code eft) for the k th mobile. 
Compute the cross-correlations px*'^ 
between eft) and the codes c(p) of all other 
p=\...k~\ users for all relative delays T. 
Store results in the look-up table. 



v 



15 



EP 1 128 564 A1 



Initiate/maintain dedicated channel with 

the hfh mobile and perform MUD by 
extracting from the table the appropriate 
correlation coefficients. 
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